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We derive a new lower bound for the ground state energy 
E^{N, S) of N fermions with total spin S in terms of binding 
energies E^{N — 1, £"±1/2) of (A'' — !) fermions. Numerical ex- 
amples are provided for some simple short-range or confining 
potentials. 

There is a persisting interest in deriving lower bounds 
on the binding energy of A'^-particle systems. Some early 
investigations were motivated by fundamental studies on 
the thermodynamic limit or the stability of matter [Q. 
Also, the search for a lower bound is rather natural once 
an upper bound is provided by variational estimates. 

Many efforts have been devoted in particular to express 
a bound on the binding energy of a iV-body system 
in terms of a {N — l)-body energy En~i with modified 
constituent mass or interaction strength. Thanks to a 
proper account for the center-of-mass motion of the sub- 
systems , the situation is now rather satisfactory in the 
zero-spin boson sector: the bound of Hall and Post |^] 
is saturated for the harmonic oscillator and approaches 
closely the exact result for many other potentials. An ex- 
tension has been derived for three and four particles with 
different masses where, again, saturation is obtained 
in the case of harmonic forces. 

The situation is far more difficult for fermions. The 
elegant bound derived by Levy-Leblond suffers from 
the fact that the energy of {N — l)-body subsystems is 
replaced by its rest energy without including its overall 
kinetic energy. As a result, saturation is never reached. 

A significant but partial improvement was obtained by 
Basdevant and Martin Q , who used subtle convexity in- 
equalities. Their bound improves that of Levy-Leblond 
in some cases, and becomes exact for the harmonic oscil- 
lator. Their approach, however, is restricted to confining 
potentials r"?, q > 1, or superpositions of such power-law 
potentials with positive weight factors. 

In this letter, we use group-theoretical techniques to 
derive a general lower bound on the A^-body energy of 
an interacting system of particles with internal degrees 
of freedom. The decomposition of the Hamiltonian is 
supplemented by considerations on the symmetry struc- 
ture of the wave function. Numerical tests are presented 



and possible generalizations are sketched, in particular 
for quantum dots, i.e., particles interacting both among 
themselves and with an external potential. 

Let us consider first N identical particles of mass m 
whose interaction does not depend on their spins. This 
corresponds to the Hamiltonian 
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where r,y = |rj — r;! is the distance between two particles. 
The case of particles with different masses is treated in 
Ref. §. 

A first decomposition of this Hamiltonian is 



N - 1 

N -2 



TO, g 



(2) 



where the constituent mass in Hj^^i is increased, since 
the kinetic energy p1/{2m), for instance, is shared by 
(N — 1) terms. The i^^ particle is absent from i/^Lj^. 
Using the variational principle with the ground state of 
the A-body system as trial wave-function leads to an 
inequality on the ground-state energies En 
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which can be rewritten differently, using the obvious iden- 
tity E{am,g) = E{m, ag)/a. This inequality is never 
saturated, because the overall translation energy of the 
(A^ — l)-body subsystems within the A^-body system is 
neglected. An improvement consists of replacing the de- 
composition by the identity 
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relating the translation-invariant Hamiltonians 
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This leads to the new inequality 

N f N ' 

EN(m,g) > I -^ m,g 
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As Nm/{N - 1) is smaller than [N - l)m/{N - 2) for 
iV > 3, the bound (H) is better than (^), since any binding 
energy in a given potential is a decreasing function of the 
constituent mass. 

By recursion, an inequality such as @ or provides 
a bound on in terms of the 2-body energy E2- In 
the case of the harmonic oscillator for bosons, the bound 
derived from ^ is saturated, whereas the bound derived 
from is by a factor -^72 smaller than the exact result 
at large N. 

Some numerical illustration of the inequality (||) for 
bosons is provided in Table |. Other examples are given 
in Ref. @. 

In the case of fermions, the bound (^) is not satisfac- 
tory. Indeed, it is not necessarily the ground state of 
which is relevant, but the lowest state with quan- 
tum numbers compatible with the symmetry of the A^- 
body state of interest. Thus looking at the structure of 
the wave function seems important to improve the in- 
equalities. 

Since our Hamiltonian acts in the orbital Hilbert space 
only, the behavior of the wave function with respect to 
the space variables has to be specified by an irreducible 
representation of the symmetry group Sn. This corre- 
sponds to a partition [A] = [Ai, A2, . . . , A„] of the N parti- 
cles defining an invariant orbital subspace spanned by or- 
thogonal states which are associated to a Young tableau 
or a Yamanouchi symbol r Both specify the trans- 
formation under permutation, and the matrix element of 
any transposition operator is known analytically. For 
example, in the case of three particles, we can have the 
following partitions 



[3], [2,1], 



[1,1,1]. 
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The first and the third ones are in the familiar sym- 
metric and antisymmetric representations of dimension 
1, respectively. The second one, of dimension 2, is 
the mixed symmetry representation spanned by the two 
states r = (2, 1, 1) or r = (1, 2, 1) which are respectively 
symmetric and antisymmetric in the exchange of parti- 
cles 1 and 2. 

Each orbital state \N, [A],r)^ must be associated with 
a similar state \N, [A], f)i for the intrinsic degrees of free- 
dom such as spin, isospin, color, etc. The coupling of 
[A] to [A] gives a symmetric representation [N] of Sn for 
bosons or a fully antisymmetric representation [1^] for 
fermions. This implies that [A],f are identical to [A],r 
when we are dealing with bosons, while for fermions, [A] , f 
correspond to a Young tableau with rows and columns 
interchanged ra]. 



The Sat Clebsch-Gordan coefficients ([A]r[A]r \[N]) and 
([A]r[A]r |[1^]) are known explicitly [^, and the iV-body 
state reads (e = 1 for bosons, —1 for fermions) 



1^, [A]), 
where 
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is the dimension of the representation [A] and a is the 
number of permutations that have to been performed 
to obtain the Yamanouchi symbol from the normal one 
@. To obtain the {N - l)-body parts of state (|), 
one has now just to remove the 7V'^ particle in each 
symbol r, thus leading to a {N — l)-Yamanouchi sym- 
bol that belongs to a different orbital symmetry [A]p = 
[Ai, . . . , Ap_i, Xp - 1, Ap+i, . . . , A„], where p is the row 
number of particle N in the original Young tableau. 
Moreover, this new symbol needs a = 1^^^ 
permutations than r to appear in the normal form. Hence 
the iV-body wave function (H) can be rewritten as 
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with p running over all the lines of [A] where a box can 
be dropped. When this is used with the decomposition 
(^ of the Hamiltonian, it leads to a new and very general 
inequality 
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If one considers, as a first application, a system of spin- 
less bosons for which the orbital Young pattern [A] is nec- 
essarily symmetric, the inequality ( pT| ) coincides with our 
previous result (|^). 

Consider now fermions with spin 1/2. The conjugate 
partition [A] of [A] is an irreducible representation of U(2) 
and contains at most 2 rows which are related to the 
spin S by the relations Ai -I- A2 = iV and Ai — A2 = 2S. 
Writing our inequality (^l]) with [A] gives the following 
result: the ground state energy of TV-spin 1/2 fermions 
with total spin S is bound in terms of binding energies 
of {N — l)-body systems with neighboring spins S" ± 1/2 
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The generalization is straightforward: the energetically 
favored state has the most symmetric space partition 
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[A] that admits a conjugate partition [A] which can be 
a representation of U(r2), where is the number of in- 
trinsic degrees of freedom. As a result, [A] has a maxi- 
mum of n rows and [A] = [fi'^ , N — z/fi] for the A^-body 
ground state with i' the integer part of N/fl. Denoting 
by -E-jv-i the binding energy of (A^ — 1) particles with 
the favored orbital symmetry [Sl^, A^ — lyU. — 1] and by 
E}^_i the energy associated to the first excited partition 
[A] = [n''-'^,n~ 1, AT - vn], the relation (|ll|) takes the 
form 



E]s{m,g) > 
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We have calculated the lower bound and the exact en- 
ergy with a selection of potentials, using the method of 
Gaussian expansion described for instance in Ref. [Q. 
The variational parameters are determined by a stochas- 
tic optimization js) . The method has proved to be pow- 
erful and reliable. 

In Table |, we show the ratio of the accurately com- 
puted energy to the lower bound for monotonic short- 
range potentials. For A^ = 3 and S = 1/2, each pair is in 
an equal- weight admixture of singlet and triplet, while 
5 = 3/2 involves only triplet states. For A^ = 4 and 
5 = (2), each 3-body subsystem has S = 1/2 (3/2). 

It appears from Table | that the bound is close to the 
exact result, especially for antisymmetric orbital wave 
functions (A^ = 3, 5 = 3/2 or A^ = 4, 5 2) and for 
deeply bound states. 

The results in Fig. [l] correspond to power-law poten- 
tials rlj with q > I. For comparison, we also display the 
bound by Levy-Leblond |l[], initially designed for large 
systems but also applicable at small A^. The decomposi- 
tion 
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expresses Hn in terms of Hamiltonians with (A^ 
dependent particles. For fermions, one gets 



A^ 
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where /jv is the cumulated energy of a system of A^ in- 
dependent fermions, a notation borrowed from Ref. [Q. 
In Eq. ([l5|), the translation energy of the (A^ — l)-body 
subsystem is neglected. 

Also shown in Fig. |l| is the bound of Basdevant and 
Martin M, who started from the identity 
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where P is the total momentum, R the center-of-mass po- 
sition and Tj- the relative kinetic energy, and succeeded to 
generalize it in the form of inequalities when the power 
q = 2 is replaced by another power q > 1. Then the 
Hamiltonian (^ or its translation-invariant part (^) with 
V cx r'^j is bounded (on both sides) by independent- 
particle Hamiltonians with a potential proportional to 



This leads to 



En > 2~ 



En<2~ N~Jn-N~E2 , 
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for 1 < g < 2. The inequalities are reversed for q > 2. 

Clearly from Fig. |l|, Basdevant and Martin always give 
the best bound near q = 2, but their result quickly de- 
teriorates when the potential departs from the harmonic 
case. Our inequality is saturated at q = 2 for A^ = 3 
fermions and for some spin configurations of Af = 4, but 
saturation is lost at larger A^. For small systems, our 
bound dramatically improves that of Levy-Leblond. 

Our result can be extended to more general Hamilto- 
nians, for instance with spin-spin interaction or with ex- 
ternal constraints. More details will be given elsewhere. 
Let us just mention that our bound can be applied to the 
"quantum dot" systems 
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once the center-of-mass contribution (a mere harmonic 
oscillator) is removed. In the case of A^ = 3 electrons 
with mass m — 1 and charge e = 1, one obtains, for 
S — 3/2 and orbital momentum and parity = 1+, a 
bound e > 0.1668 if the oscillator frequency is = 0.01 
and e > 54.968 ii uj — 10, to be compared to the exact 
values e = 0.1680 and e = 54.973. For A^ = 4, 5 = 2, 
= 0", one obtains e > 0.297 if w = 0.01 and e > 
84.895 if = 10, close to the exact values e — 0.299 
and e = 84.907. The bound is better for large uj, as the 
system becomes a pure oscillator. 

To conclude, we have derived a lower bound on the 
energy of A'^-fermion systems, which is independent of 
the shape of the potential. For small systems, it improves 
significantly a previous bound expressed in terms of ( A^ — 
1) independent fermions. 

Saturation is not always obtained for the harmonic os- 
cillator. When our bound is iterated to express the A^- 
body energy in terms of the 2-body energy, the correct 
large A^ behavior is not reached. Thus there is still much 
room for improvement. 
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TABLE L Results for Yukawa (Y), Gaussian (G) and expo- 
nential (E) potentials. The range parameter is set to unity by 
rescaling. The quantity shown is the ratio of the computed 
energy to the bound. In the first line, the entries are: the 
number of particles, the boson (B) or fermion (F) character, 
and the total spin. The 4-body energy is compared to the 
3-body one, as per Eq. (^) for bosons and (^^ for fermions. 
The strength g refers to the Hamiltonian for A'^ = 3 particles. 



iV = 3, 5* = 1/2 
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FIG. 1. Comparison of various lower bounds for A*' = 3 or 4 
fermions with total spin S , in the case of power-law potentials 
r'. The quantity shown is the ratio of the exact energy to our 
bound (thick line), to that of Basdevant and Martin (dotted 
line) and that of Levy-Leblond (thin line) . 
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